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I start at CERN in 2008. The experiment is primarily dedicated to B-Physics and hence requires 
' precise vertex reconstruction. The sihcon vertex locator (VELO) has a single hit precision of better 
. than 10 /xm and is used both off-line and in the trigger. These requirements place strict constraints 
on its alignment. Additional challenges for the alignment arise from the detector being retracted 
' between each fill of the LHC and from its unique circular disc r/<j) strip geometry. This paper 
■ describes the track based software alignment procedure developed for the VELO. The procedure 
. 5^ I is primarily based on a non-iterative method using a matrix inversion technique. The procedure is 
demonstrated with simulated events to be fast, robust and to achieve a suitable alignment precision. 



Alignment procedure 
of the LHCb Vertex Detector 

S. Viret, C. Parkes, M. Gersabeck 

Department of Physics and Astronomy, University of Glasgow 
University Avenue, Glasgow, G12 8QQ, United Kingdom 

Abstract 

LHCb is one of the four main experiments of the Large Hadron Cohider (LHC) project, which wiU 
start at CERN in 2008. The experiment is primarily dedicated to B-Physics and hence requires 
precise vertex reconstruction. The silicon vertex locator (VELO) has a single hit precision of better 
than 10 fj,m and is used both off-line and in the trigger. These requirements place strict constraints 
on its alignment. Additional challenges for the alignment arise from the detector being retracted 
between each fill of the LHC and from its unique circular disc r/(j) strip geometry. This paper 
describes the track based software alignment procedure developed for the VELO. The procedure 
is primarily based on a non-iterative method using a matrix inversion technique. The procedure is 
demonstrated with simulated events to be fast, robust and to achieve a suitable alignment precision. 
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1 Introduction 

LHCb is the dedicated heavy flavour physics experiment at the LHC. The physics goals are 
critically dependent on the performance of the precision vertex locator (VELO). Whilst the 
intrinsic VELO sensor single hit resolution is 5 — 10 fim for all tracks in the acceptance, the 
sensors also have to be retracted by 3 cm from the unstable LHC beams while the machine is 
fllled. As a result of these circumstances and the unique geometry of the detector, described in 
Section 2, the VELO has particularly demanding alignment requirements. 

The VELO detector has been assembled with a high precision and a detailed metrology has been 
performed. However, the only possible method for correcting misalignments during data taking 
is through a track-based software alignment procedure. A new alignment will be performed 
after each re- insertion of the VELO, in order to check the alignment and ensure optimal data 
taking. 

The alignment procedure is based on techniques which are widely used in particle physics, a 
comprehensive review of particle physics alignment methods is available in Ref. ([ll). Two types 
of techniques are applied here, one is based on matrix inversion and the other on the fltting of 
functional forms to residual distributions. The LHCb VELO alignment is however unique due 
to the module and detector geometry and the frequent mechanical retraction of the detector. 

The VELO alignment algorithm is divided into three phases: the relative alignment of the 
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sensors in each module; the relative alignment of the modules in each half of the detector; and 
the relative alignment of the detector-halves. The algorithms are described in Section 3. The 
results obtained using simulated events are described in Section 4. A summary and conclusions 
are given in Section 5. 

The performance of the first two phases of the alignment procedure were also tested on data 
during a VELO beam test. The results of this study are presented in Ref. (Q) and verify the 
simulation results presented here. 



2 The VELO Alignment Context 



The VELO consists of two retractable halves, each with 21 modules, as shown in Fig. 1. Each 
module contains two semi-circular silicon strip sensors that measure r (radial) and (azimuthal 
angle) co-ordinates, respectively. The strips on the measuring sensor have a stereo angle off- 
set. A hole up to a radius of 7 mm in the centre of the sensors allows the two proton beams to 
pass through the VELO, with their interaction point being close to the start of the detector. 




Fig. 1. Layout of the modules in the LHCb vertex detector (left). The right figure describes, in the 
VELO reference frame, the rotational degrees of freedom used throughout this paper. Rotations around 
the different axis are applied in the following order: Z (A^), Y' (A^), and X" (A^). 

The VELO sensors are operated in vacuum and separated from the primary beam vacuum of 
the LHC by a thin (250 /xm) aluminium foil. The two halves of the VELO will be retracted by 
3 cm during each LHC fill until stable beam conditions have been established. Further details 
of the VELO detector design can be found in Ref. (j^l). 

No magnetic field is applied in the VELO active area, and the residual effects of the LHCb 
dipole magnet are small. Thus the tracks reconstructed in the detector can be represented by 
straight lines in the VELO at first order. The effect on the alignment of using straight lines, 
and hence neglecting the magnetic field and multiple scattering effects, in the VELO have 
been studied in simulation. Magnetic field effects could affect the alignment result only when 
a low momentum track sample with a large charge asymmetry is used. Those effects may be 
neglected when they are minimized by using an equal number of positively and negatively 
charged particles. Furthermore, one can use the LHCb Trigger Tracker (0) information (also in 
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an un-aligned environment) to select a track sample with equal numbers for both charges. 

The successful operation of the LHCb trigger puts tight constraints on the VELO assembly 
precision. The LHCb trigger system relies on vertexing but, for speed reasons, the VELO 
pattern recognition algorithm has to be performed initially in the r-z projection. This requires 
that the strips on the R sensors accurately describe circles around the beam position. 

Moreover, remaining misalignments that are not corrected for in software could have significant 
effects on the trigger performance. For example, a rotation of 0.5 mrad around the Y-axis of 
one detector can reduce the trigger efficiency by 30 % for the Bg K^K~ channel (jsj). 



3 The VELO Alignment Procedure 



The alignment of the LHCb VELO proceeds in a number of stages: precision assembly; metrol- 
ogy; and the software alignment described in this paper. 

A survey of the VELO system has been performed. The relative position of the R and $ sensors 
on a single module have been measured to an accuracy better than 5 //m, and no significant 
curvature of the sensors was observed. The relative position of the VELO modules on the 
VELO half bases has also been determined with a similar precision. The survey is of particular 
importance for determining degrees of freedom that can be constrained less well from tracks, 
such as the relative z-positions of the VELO modules. Moreover, the survey values have been 
taken as the initial VELO alignment conditions: the software alignment will be seeded with 
these values. 

The VELO software alignment procedure is divided into three phases. These phases are moti- 
vated by and refiect the mechanical construction of the system and the different time periods 
on which re-alignment is anticipated. 

• The first phase performs a relative alignment of R and $ sensors within one VELO module (@). 
The R and $ sensors are glued together onto the same hybrid, thus making their relative 
alignment highly stable. It is anticipated that the alignment of this structure will not need 
to be repeated as frequently as for the other phases. The alignment technique is based on an 
iterative fit of the distribution of track residuals across the surface of the silicon sensor. 

• The second phase is an internal alignment of the modules within each VELO detector-half (0). 
The two VELO detector halves are mechanically separate and can be retracted, partially 
inserted or fully inserted. The VELO halves will be retracted each fill and hence, at least 
as a cross-check for relative module movement, the alignment will be performed each fill. 
The alignment technique is based on track residuals applying a non-iterative method using 
a matrix inversion technique. 

• The third phase is the relative alignment of the two halves with respect to each other (jsj). 
This will be performed each fill to cross-check and improve upon information provided from 
the mechanical insertion system. The alignment technique uses the same technique as the 
module alignment but requiring tracks that pass through both VELO halves and on a similar 
technique that is applied with vertex constraints. 

These three alignment algorithm phases are described in more detail in the following sections. 
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3.1 Relative Sensor Alignment 



A VELO module contains an R and a $ sensor glued onto the same hybrid, so that the 
sensors are back-to-back. The first step of the alignment procedure is to determine the relative 
alignment of the R and $ sensors. This relative alignment is already known to high accuracy 
from a mechanical survey but can be further improved by applying the following procedure to 
determine the critical relative x and y translations of the sensors. 

As the R and $ sensors on each module have to be treated separately, the residuals with respect 
to strip hits have to be considered. The unique VELO R/^ sensor geometry means that the 
hnearized non-iterative approach applied to space-points in the subsequent alignment phases is 
not applicable. Instead, an iterative alignment procedure that fits the distribution of residuals 
as a function of the sensor azimuthal angle has been developed. This procedure is partly inspired 
by the method developed for the SLD vertex detector (Q). 

The observed signals on sensor strips are used to determine the best estimate of the hit position. 
The hits on the sensors (other than those in the module under study) are fitted to produce 
tracks. These tracks are extrapolated to obtain the track intercept point in the sensor under 
study. The unbiased residual is then defined as the perpendicular distance between the track 
intercept point and the line parallel to the strip at the observed hit position. Consequently, 
these residuals are sensitive to sensor misalignments perpendicular to the strip direction. Since 
the direction of the strips changes as a function of the azimuthal angle the residuals thus have 
a well defined sensitivity to translational misalignments in x and y (see Fig. 2). 
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Fig. 2. Influence of misalignments on residuals of R and $ sensors. The misalignment shown is a 
translation both along the negative x and y direction. 

It can easily be shown that the relation between residuals {tR/q,) and misalignments in the 
sensor plane (Ax, Ay, A^) is given by 



Ci? = -Ax cos (j)track +Ay siu ^i^acfc 



(i? sensor) 



e$ = Ax sin + cos (ptrack +^-yrtrack (* seusor) 



(1) 



where (ptrack is the azimuthal angle of the extrapolated track position. A^ describes a mis- 
alignment in the form of a rotation around the z axis, which translates into a shift in by 
multiplication with the radial co-ordinate of the extrapolated track in the sensor plane. As this 
does not contain any dependence it is sufficient to leave it as a free parameter in the fit. 
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In practice the alignment constants are determined by an iterative fitting algorithm. For each 
iteration, the unbiased residuals for the R and $ sensor are plotted against the azimuthal angle. 
The misalignments are determined through a fit to a binned distribution of the mean residual 
values as function of (p. 

An additional complication arises from the fact that the strips on the VELO $ sensors are not 
exact radial lines. The $ sensors are divided into an inner and an outer region within which 
the strips are tilted through a stereo angle of 20° (inner) and —10.35° (outer). This can be 
compensated for by making the replacement in Eq. 1 

where f3 is the stereo angle and 0„tm is the coordinate at the minimal radius of the strip. In 
addition to preserving the simple relation between the residuals and the alignment constants, 
this transformation also allows the fit to be made to data from the inner and outer region of 
the $ sensor simultaneously. 

This method is primarily used to determine the relative x and y translations of the sensors 
in a module. However, to improve the fit convergence, the value for A^, the rotation around 
the z-axis of the $ sensor, is also determined. This is done using an equivalent approach as 
above, however now exploiting the distribution of residuals as a function of the track radial 
co-ordinate. The value for is then determined as the slope of fit with a linear function. 
However, this relative z rotation of the modules will be extracted to a higher precision using 
the technique described below in step two of the alignment procedure. 

The fitted tracks used in the procedure rely upon the current estimate of the alignment constants 
of the sensors. Hence it is necessary to iterate this procedure, updating the fitted tracks as the 
alignment constants are improved. 

As the unbiased residuals are determined from a track fitted while excluding hits on both sensors 

under study, the influence of module to module misalignments vanishes to first order when 
determining the relative translational misalignment of the sensors. The only misalignments 
that cannot be determined with this method are common sensor to sensor misalignments of all 
modules, and sensor to sensor misalignments that have a linear dependency on z. 

The results of applying this first step of the alignment procedure to simulation events are given 
in section 4.1. 

3.2 Relative Module Alignment 

The second step of the alignment procedure is to perform the relative alignment of the 21 
VELO modules within one detector-half. This is based on a non-iterative method using a 
matrix inversion technique to minimize a function. 

The is produced from the residuals between the tracks and the measured clusters. The mea- 
sured residuals arc expressed as a linear combination of the track parameters and the alignment 
constants. Each straight line track has four parameters (two slopes and two intercepts) denoted 
niocai- The total number of translational and rotational degrees of freedom of all of the modules 
is Ugiobai alignment constants. The track fitting and alignment problem can then be expressed 
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as a system of ritotai equations, where ritotai is given by 

^total ^local ' ^tracks ~l~ ^global- 



Clearly, the size of the system scales with the number of tracks used in the alignment. As 
discussed later, about 20,000 tracks are necessary to obtain a detector-half alignment of the 
required accuracy. Hence, this implies solving a system of over 80,000 equations, which is a 
computationally challenging task. However, this problem can be reduced to the size of Ugiobai 
(around 100 for the VELO) using the technique of matrix inversion by partition, performed 
by the program Millepedqlj (lO). This algorithm has already been successfully used in many 



particle physics experiments (see for example (llll: Il2l)) 



As stated previously, VELO /2-sensor strips are semi-circular, thus making i?-sensor hit infor- 
mation non-linear. However, in order to use the Millepede approach it is mandatory to establish 
a linear relation between the hit residuals and the misalignments. Here, this is obtained by pro- 
ducing space-points at the VELO module level as described below. This is the natural choice 
for the VELO as a module containing an R and a $ sensor is a independent mechanical object 
whose internal alignment is expected to be relatively stable. 

A particle passing through a VELO module (which is assumed to have its R and $ sensors al- 
ready internally aligned) gives a non-linear set of coordinates: (r, z(i?sensor)) and {(p, sensor))- 
This system is transformed into an (x, y, z) space-point by projecting the information onto 
the R sensor. As the co- ordinate of the particle can change slightly in the gap between 
the R sensor and the $ sensor, a correction based on the currently estimated track slopes is 
applied. Since the correction is small, and the track slope estimates are not significantly biased 
by the expected module misalignments, the non-linear effects introduced by this procedure are 
negligible. 

An (r, (f)) cluster pair is thus transformed into an (x, y, z) space-point as follows: 

X = r • COs(0corr) 

y = r • sin(0corr) 

2 '^{Rsensor) 



with 



~l~ (ptrack(^Rsensor) 'Ptracki^^ sensor) • 



The linear relation between the residuals and the alignment parameters for the VELO modules, 
which is derived in Ref. (0), is given by: 

= -Ax +y +a-(A, +x- A;3 + y A„) 
ey = -Ay -X • A^ +c • (A^ + X ■ A/3 + y • A„) 



where x, y, z represents the measurement values, and a, c are the slopes in the XZ, YZ planes 
respectively of the track considered. The Millepede technique is then used to simultaneously 

^ a C++ translation of the FORTRAN Millepede program has been implemented by one of the authors 
of this paper (13) • 
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extract the alignment constants (Ax, Ay, A^, A^, A/3 and A^) of each of the 21 VELO 
modules. 

An important component of the alignment procedure is the selection of the track sample used. 
In order to ensure an optimal population of the final matrix, a mixture of tracks coming from 
the primary interaction point and a complementary set of tracks from the beam halo or beam- 
gas interactions will be used. A specific pattern recognition algorithm has been developed in 
order to select these events (flih . 

The 'weak-modes', i.e. deformations which are difficult — if not impossible — to unfold with 
tracks, are extensively constrained. Two categories of 'weak-modes' have been considered. 

( TRANSLATIOhTj ( ROTATION ] 




[ SHEARING J [ SCALING ] 



Fig. 3. Linear deformations impossible to constrain via VELO module alignment. 

Firstly, track residuals are insensitive to linear transformations of the whole detector-half {e.g. 
translation along an axis, or rotation around an axis). The different possible effects are illus- 
trated pictorially in Fig. 3. In three dimensions, this leads to 12 possible global deformations. 
However, due to the structure of the VELO, we can neglect three of them at first order: shearing 
in the XY plane, and scaling of the X and Y axes. In addition, due to the VELO geometry, 
it will be difficult to distinguish XZ and YZ shearings from Y and X rotations. We choose to 
constrain only the XZ and YZ shearings, as they are easier to add into our linear system. 

Hence we constrain 7 possible deformations: Z axis rotation, X, Y, Z translations, XZ and YZ 
shearings, and Z axis scaling. These degrees of freedom are fixed in order to prevent these 
movements occurring during the module alignment. Within the Millepede framework, this is 
achieved through Lagrange multipliers, i.e. through the addition of constraint equations to the 
global system. 

The second class of 'weak modes' are those that do affect the track residuals but to which there 
is greatly reduced sensitivity. This is, for example, the case for module rotations around the 
X-axis, or around the Z-axis with an angle proportional to Z. Even with the best possible track 
sample the effect on the of these modes is comparatively weak, hence they are also weakly 
constrained by any alignment procedure. Using the survey information as a starting point 
for the alignment will provide a strong constraint against those special modes. In addition 
to that, and in order to avoid non-physical variations during the alignment procedure, each 
alignment constant is controlled: 'penalty-terms' are added to the to minimize. These extra- 
terms prevent the variation of the alignment parameters (w.r.t. the initial survey value) to be 
significantly larger than the metrology precision. 
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The results of applying this second step of the alignment procedure to simulation events are 
given in section 4.2. 



3.3 Detector-halves Alignment 

The final step of the alignment procedure is to perform the relative alignment of the two 
detector-halves. Clearly, the module space-point residuals within a detector-half are not sensi- 
tive to the misalignment of the whole half. Thus, observables connecting the two detector-halves 
have to be used. 

Tracks that pass through both VELO detector-halves, referred to here as overlap tracks, provide 
a powerful constraint for the relative positioning of the two detector-halves. These tracks are 
required to have at least one space-point {i.e. an (r, 0) cluster pair) in each VELO detector- 
half. The VELO was specifically designed to obtain such a class of tracks: the modules from the 
VELO left and right-halves are offset in z and interlace by a small distance when closed. The 
overlapping area between the two VELO detector- halves, when fully inserted, corresponds to 
2% of the active surface area of the sensors. 

Having selected these overlap tracks, the alignment proceeds by residual minimization using the 
same technique as the module alignment. An equivalent of Eqn. 2 is thus required in terms of 
the relative alignment constants of the two detector-halves. Assuming that the alignment of the 
modules in each detector-half has been already corrected, then the equation for the residuals 
becomes (see Ref. (js])): 

fe. = -A^ + z^Af + yA^ + a.(Af + x-Af + yAf) 
[ey = -Af + z^Af-x•A^ + c•(Af +x-A^ + yAf) ' 

where A^, Ay , Af , Af , A^, A;^ are the detector-half degrees of freedom (position of one 
half with respect to the other), and z^ the z position in the detector-half frame of the module 
in which the space-point hit was recorded. 

The matrix inversion technique then allows these six relative alignment constants of the detector- 
halves to be determined, assuming the VELO is fully inserted for physics data taking. 

However, during insertion and during the commissioning phase an alignment may be required 
with the VELO in the retracted position. In this case the rate of overlap tracks becomes 
very small and an alternative technique is required. The detector-half alignment can then be 
performed using vertices as constraints. Tracks from the primary interaction point can be used if 
available or beam-gas interactions occurring in the VELO vacuum tank. By fitting for the vertex 
separately inside each of the two detector-halves, one can obtain the misalignment between the 
two halves. 

Assuming that one box is fixed {i.e. using it to define the co- ordinate system) the reference 
vertex position can be determined from this half. The vertex position found using the tracks 
passing through the other detector-half is related to the misalignment between the two halves 
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via the following relations: 

= k + a,- + - tti • Af - ■ Af 
^v^ = d^ + Q- + Af - c, • Af - ■ Af 

where a^, bi, Ci, di are the parameters of the i^^ track used to fit the vertex, (ff , ff , ff ) the vertex 
position in the 'moving' half, and Af denote the detector-half degrees of freedom (position of 
one half with respect to the other). It is easy to relate this formulation to the residual formalism 
used previously. Indeed, the residual in this case will be the differences between the vertices 
positions fitted within the two halves: 



where {v'^^ ^v"^^^ ^v^^^) defines the vertex position in the fixed half, i.e. the reference value. 

Having determined these relations, the same technique can again be used to determine the 
alignment constants: even though we are now using vertices rather than track residuals the 
Millepede framework is still applicable as it does not depend on the object that is being fitted. 

More details on this novel technique of alignment with vertices can be found in Ref. (jsl). The 
vertex fitting method can also be used to determine the VELO detector-halves' positions with 
respect to the beam, when using vertices from the primary interaction point. 

The results of applying this final step of the alignment procedure to simulation events are given 
in section 4.3. 



4 Simulation Studies 

A simulation of 200 samples of 25,000 events each has been produced and propagated through 
the LHCb software. Each sample, which comprises a mixture of 5,000 minimum bias events 
[k, 100,000 tracks from primary vertex interactions) plus 20,000 beam-halo like events, was 
produced with a different set of alignment constants. The misalignment values have been ran- 
domly chosen within a Gaussian distribution centered on zero and with resolutions based on 
construction and survey accuracies (defined in Ref. (j^)). All the module and detector-half de- 
grees of freedom have been misaligned. The initial step, the relative sensor alignment, was 
tested using different event samples which are described in the first part of this section. 

^.1 Relative Sensor Alignment Results 

The sensor alignment method has been tested with 10 samples of randomly generated misalign- 
ments. All sensors have been misaligned individually, thus generating a scenario equivalent to 
simultaneous module to module and sensor to sensor misalignments. Each of the 10 samples 
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consists of 20,000 tracks with small slopes, thus passing through all sensors of one VELO-half 
and evenly distributed across the sensor surface. Typically three iterations of the alignment 
procedure are required to obtain the best resolution. 

Fig. 4 shows the generated and the remaining misalignments after all iterations. The resolution 
on the relative x and y translation of the sensors of one module is 1.3 fim, i.e. a significant 
improvement over the survey precision. The performance of this algorithm has also been demon- 
strated with beam test data and is reported in Ref. 
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Fig. 4. Misalignment values before (□), 
and after (■) sensor relative alignment. 



Fig. 5. Misalignment values in one de- 
tector-half before (□), and after (■) in- 
ternal alignment. 



4.2 Relative Module Alignment Results 



The internal alignment of the modules in each detector-half is primarily sensitive to translations 
of the modules in the X and Y directions and rotations around the Z-axis. Results on the 
alignment of all modules in a detector-half, for one particular sample, are presented in Fig. 5. 
In Fig. 6 the alignment constants for 200 event samples, before and after correction, are shown. 
Resolutions on the X and Y translation alignment parameters of 1.1 fim and on rotations 
around the 2;-axis of 0.12 mrad are obtained. About 20,000 tracks per detector-half were needed 
to obtain this accuracy. The alignment resolutions are well below the intrinsic detector hit 
resolution. The performance of this algorithm has also been demonstrated with beam test data 
and is reported in Ref. (j^. 

Concerning the 'non-linear' degrees of freedom, the observed sensitivity is as expected worse 
than for the other parameters. However some results were obtained for the modules which are 
close to the interaction region, i.e. where track slopes are larger. Restricting the study to these 
stations (1 to 14), one obtains a reasonable sensitivity to (28 /im) and a fair sensitivity to 
Aa and (0.8 mrad and 1.1 mrad respectively). This sensitivity is worse than the survey 
precision, but will provide a cross-check of this survey information. 

As this algorithm is in general run independently of the relative sensor alignment algorithms 
its performance has been evaluated separately. In the presence of relative misalignments of the 
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sensors on a_given module the module's position will be aligned to the average position of the 
two sensoro- 



4.3 Detector-halves Alignment Results 



Although the three alignment steps can be performed independently, in practice it is expected 
that steps two and three will be run consecutively. Hence, the results presented in this section 
are for the realistic case of performing both of these alignment steps on misaligned samples. 
The tracks are refitted after the module alignment procedure in order to update the track 
parameters. The results presented here have been obtained with about 300 overlap tracks. 

The results of the study are presented on Fig. 7. The resolution on the X and Y translation 
alignment parameters is 12 /im for x and y translations, and the resolution on the x and y tilts 
is 36 /irad. 
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Fig. 6. Misalignment values before (□), 
and after (■) internal alignment, for all 
configurations and all stations. 
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Fig. 7. Misalignment values before (□), 
and after (■) detector halves align- 
ment, for all configurations. 



As in the case of the module alignment, some of the degrees of freedom are more difficult 
to constrain. In the detector-half alignment these weakly constrained motions are the ones 
related to the Z-axis: rotation around and translation along. The relative rotation around Z 
between the two detector-halves is constrained using the overlap tracks. Translations along 
Z are estimated through the vertex fitting technique. The vertices are fitted separately from 
tracks in each detector-half and the misalignment determined by comparing the Z positions, 
from this a 40 fim resolution was obtained. 



This requires a track sample with a sufficiently flat distribution in 
used for VELO alignment. 



which is given for the samples 
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5 Conclusion 



A software alignment method for the LHCb vertex locator has been developed. This procedure 
performs the alignment in three steps, all using track residuals. First, a relative alignment of 
the R and $ sensors within each VELO module is performed. This is followed by an internal 
alignment of the VELO modules within each detector- half. Thirdly, since the detector-halves 
are moved between each LHC fill, a final step is required in order to align the detector-halves 
with respect to each other, thus providing a fully internally aligned VELO. 

Due to the VELO module design {R and $ sensors are glued onto the same hybrid and precisely 
surveyed), it is foreseen to perform the first stage with a much lower frequency than the two 
other ones and to use it for the annual data reprocessing. The alignment of the relative sensor 
position achieves a precision of 1.3 fim for x and y translations. 

The final two stages of VELO alignment strategy use the Millepede program, which enables the 
alignment to be performed in only one pass; this is to be compared with classic minimization 
methods which require many iterations to provide their result. This technique allows the pro- 
cessing of module and VELO-half alignment within a few minutes on a sing le CP\M, assummg 
that an appropriate data sample is available. 

For the internal alignment of the modules in a detector-half a 1.1 /im precision has been 
achieved on the relevant translational degrees of freedom {i.e. along X and Y axes), and a 
0.1 mrad accuracy on the rotation around the Z-axis. These values are well below the expected 
detector intrinsic resolution. Future work will be undertaken to enable an efficient selection of 
the required tracks in the LHCb trigger. 

The results of the detector-half alignment procedure show that the tracks that pass through the 
overlap region between the two VELO detector-halves provide a very strong constraint. With 
only a few hundred tracks accuracies of 12 /xm for x and y translations, and 36 /xrad for X and 
Y tilts, are obtained. These results are well within the system requirements. 

This paper has described the algorithms developed to perform the alignment of the LHCb 
VELO and demonstrated with simulation samples that this approach provides the performance 
required so that misalignment effects will not adversely affect the LHCb physics programme. 
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